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1.  Introduction  and  Summary 

Inequalities  play  an  important  role  in  reliability  theory.  One  example 
is  the  following:  Let  F be  the  lifetime  distribution  of  a component  or 

system  with  failure  rate  function  r(*)  . Suppose  inf  r(t)  = rQ  > 0 and 

sup  r(t)  = r^  < 00  . If  Fq  and  F^  are  cdf's  of  exponential  random  vari- 
ables with  means  r^  and  r^  respectively,  then  Fg(x)  - F(x)  - F^  (x) 
for  x < 0 . Such  inequalities  are  valuable  in  situations  where  the  lifetime 
cdf  F is  unknown  but  it  is  possible  to  know  or  establish  bounds  on  the 
failure  rate  function. 

If  a component  is  instantaneously  replaced  or  repaired  when  it  fails, 
the  sequence  of  failure  times  forms  a renewal  process:  Let  X.,  i=l,2,... 

be  a sequence  of  i.i.d.  lifetime  random  variables.  Let  Sp  = 0 , S = X , 

S2  = X1  + X2  ' S3  = X1  + X2  + X3  ' "■  ' then  * Si'  i=cM<2>---  } is  a 
renewal  process.  The  purpose  of  this  note  is  to  present  bounds  for  the 

renewal  process  S similar  to  the  example  of  the  preceding  paragraph, 

namely:  Suppose  that  the  distribution  of  X has  failure  rate  function  r(*) 

such  that  inf  r(t)  = r > 0 and  sup  r(t)  = r < 00  . We  show  (Corollaries 
t 0 t 1 

1 and  2)  that  it  is  possible  to  construct  Poisson  processes  { T*?,  i=0,l,2, 
...  } and  { T^,  i=0,l,2,...  } with  rates  rQ  and  r respectively,  such 
that 


{ T°  T°  T° 

1 0'  1'  2' 


} " { V V V 


1 1 T0'  1'  T2'  •••  J 


(1.1) 


with  probability  1 , i.e.,  the  failure  times  of  the  respective  processes 
are  subsets  of  one  another.  (We  shall  actually  construct  nonhomogeneous 
Poisson  processes  which  satisfy  (1.1)). 
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It  is  possible  to  use  this  to  compute  bounds  on  expected  maintenance 
costs,  for  example.  This  result  is  probably  also  of  general  interest  as  a 
contribution  to  the  theory  of  Poisson  processes  and  renewal  theory  and  thus 
may  be  useful  for  establishing  bounds  in  queueinq  theory  or  inventory  theory. 

It  is  also  possible  to  obtain  a comparison  result  (Theorem  3)  for 
certain  alternating  renewal  processes  which  arise  when  repairs  or  replace- 
ment are  not  instantaneous.  Barlow  and  Proschan  (1976)  consider  such 
processes.  Theorem  3 and  its  proof  give  a slightly  different  and  more 
general  approach  to  their  problem.  (We  assume  NWU  repair  times;  they  assume 
DFR. ) 

There  is  a growing  literature  on  comparison  results  for  stochastic 
processes.  There  are  two  general  approaches:  stochastic  inequalities 

[4,8,11,15,24]  and  almost  sure  inequalities  [13,14,22].  (The  list  of 
references  is  by  no  means  exhaustive.)  Pledger  and  Proschan  (1973)  , Ross 
(1974) , Keilson  (1974) , and  Barlow  and  Proschan  (1976)  have  applied 
stochastic  comparison  techniques  to  reliability  of  maintained  systems. 

The  comparisons  derived  in  this  paper  will  be  of  the  "almost  sure"  variety. 
These  results  are  special  for  renewal  and  Poisson  processes  and  are  derived 
from  first  principles,  not  from  any  theories  developed  in  the  above-cited 
literature,  for  example  [l3,22]  do  not  apply.  The  above  theories  tend  to 
concentrate  on  stochastically  monotone  Markov  processes  [4].  We  take  more 
of  a point-process  point  of  view  in  this  paper,  focusing  on  sample  paths 
and  intensities  rather  than  transition  probabilities.  Our  approach  is 
constructive  and  indicates  some  interesting  relationships  involving  the 
failure  rate  functions;  it  may  have  pedagogical  value. 

We  shall  use  the  following  relationship  between  stochastic  and  almost 
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sure  inequalities  for  the  univariate  case:  Let  X and  Y be  real-valued 

random  variables  with  cdf's  F and  G,  then  X is  stochastically  less 

St 

than  Y (denoted  X Y ) if  F(z)  > G(z)  for  all  z.  Let  U and  V 

be  real-valued  random  variables  defined  on  the  same  probability  space 

as 

(fi..F,P),  then  U is  almost  surely  less  than  V (denoted  U < V ) if 
P ( Cl)  e ft  : U (w)  V(u)  ) = 1 . Let  X = U moan  that  the  random  variables 
X and  U have  the  same  distribution. 


Lemma  1.  Let  X and  Y be  real-valued  random  variables.  If  X $ Y, 

then  there  exist  a probability  space  (i2,F,p)  and  random  variables  L1 

D D as 

and  V defined  on  it  such  that  U ■■  X,  V = Y , and  U < V. 


Proof : Let  be  [c,lj,  f be  the  Borel  sets  and  P be  Lebesgue 

measure.  If  X and  Y have  cdf's  F and  G , let  U = F ^ and 
V = G_1. 


(This  result  can  be  extended  to  random  vectors  and  random  functions  [ 1 2 ] . ) 


2.  Poisson  Bounds  on  Renewal  Processes 

In  this  section  the  existence  of  (dependent)  nonhomogeneous  Poisson 

and  renewal  processes  satisfying  (1.1)  on  a common  probability  space  is 

verified.  We  will  use  two  representations  of  point  processes  on  [o,00)  : 

the  counting  process  N and  the  partial  sum  process  S . N(t)  = # points 

in  (0,t].  S = inf  { t : N(t)  < n } , n = 1,2 There  will  be  a 

n 

point  at  0 : = 0 ; this  is  the  zeroth  point  and  is  not  counted  by  N - 

The  necessary  background  facts  about  homogeneous  Poisson  processes,  renewal 


processes,  failure  rates  and  other  topics  which  will  be  used  without 


-4- 


reference  can  be  found  in  most  introductory  stochastic  process  texts 
[l,6,9,19,22].  The  necessary  facts  about  nonhontogeneous  Poisson  processes 
on  [o,00)  are  presented  below. 


Definition  1 Let  r(*)  be  a real-valued  function  on  [o,00)  which  is 
integrable  over  bounded  intervals,  then  { N(t),t  S 0 } is  the  counting 
process  representation  of  a nonhomogeneous  Poisson  process  with  intensity 
function  r if 

i)  N (0)  = 0 

ii)  N has  independent  increments 

iii)  for  0 < s < t , N(t)  - N(s)  is  a Poisson  random  variable  with 
,t 

E(  N (t)  - N(s)  ) = / r(u)du  . 

s 

ft 

We  define  A(t)  = I r(u)du  as  the  cumulative  mean  function  of  the 
Poisson  process  with  intensity  r(*)  . 

Lemma  2.  Let  { N(t),  t S 0 } be  a nonhomogeneous  Poisson  process  with 

intensity  function  r(*)  . Let  T = inf  { t : N(t)  = n } , n = 0,1, 

n 

2,...  . Given  that  N(t)  = k , then  T^'  T2'  •••  • T}c  are  jointly  dis- 
tributed as  the  order  statistics  of  a sample  of  size  k from  the  cdf 
A(s)/A(t)  , OSsSt  . The  density  of  this  conditional  joint  distribution  is 


ft(  tl'  t2' 


' ^1. 


k ) = k! 


r (t . ) 
1 


i=!  ATtT 


0StlSt2S 


< tk  < t 


Proof : Note  that  the  conditional  distribution  of  the  interarrival  time 

T.,.  - T.  given  T.  = t.  is 
l+l  i 11 
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p{  T , - T.  > I T.  = t.  } 
l+l  l 1 1 l 

= P{  N(ti+X)  - N(tJ  = 0 | N(tJ  = i } 

t.+x 

= exp(  - / r(s)ds  ) 
t. 

l 

Thus  a version  of  the  density  is 

t.+x 

fT  -T  It  (xlti)  = r(ti+*>  exP(  " J r(s)ds  ) . 
i+1  i i ti 

The  conditional  density  f ( t^ , . ..,  t^  | k ) can  be  computed  now 

exactly  as  is  done  in  the  proof  of  the  analogous  result  for  homogeneous 
Poisson  processes,  see  [ 19]  p.17,  [l]  p.67,  or  [9]  p.126. 

Lemma  3.  Let  r(*)  be  right  continuous,  then  {N(t),t^o}  is  a nonhomo 
geneous  Poisson  process  with  intensity  function  r(*)  , if  and  only  if, 
for  t > 0 , 


h'1  P{  N (t+h)  - N(t)  =1  | N (s) , s < t } = r(t)  (2 


lln  h_1  N ( t+h)  - N(t)  > 2 j N(s) , s < t } = 0 


almost  surely,  and  N(0)  =0  . 


Proof : If  N is  a nonhomogeneous  Poisson  process  N(t+h)  - N(t)  is 

t+h 

independent  of  { N(s),  s ^ t } and  is  Poisson  with  mean  r(u)du  . 

t+h  t+h 

P{  N(t+h)  - N(t)  =1  }=exp(  - /t  r(u)du  ) f r(u)du  = hr(t)  + o (h)  . 
This  verifies  (2.1)  ; (2.2)  is  verified  similarly. 
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Convcrsely,  (2.1)  and  (2.2)  imply  N is  nonhomoqeneous  Poisson.  Part 
(iii)  of  Definition  1 can  be  verified  by  solving  the  Kolmogorov  forward 
differential  equations  as  in  the  homogeneous  case,  [9]  p.24,  [20]  p.118. 

Part  (ii)  of  Definition  1 follows  because  the  limits  are  independent  of 
the  entire  past  history,  { N(s),  s < t } : this  implies  the  Markov 

property  which  in  turn  implies  independent  increments  in  this  case. 

Note  that  if  (2.1)  and  (2.2)  were  weakened  by  conditioning  over 
{ N(t)  } instead  of  { N(s),  s S t } then  (unless  the  Markov  property  is 
assumed)  N is  not  necessarily  a Poisson  process.  Similarly  part  (ii)  of 
Definition  1 cannot  be  weakened  (to  pair-wide  independence,  for  example). 
Shepp  (see  [5]) and  Szasz  (1970)  present  counterexamples  for  the  homogeneous 
case;  see  Renyi  (1967)  also. 

We  can  now  present  the  first  comparison  theorem. 

Theorem  1.  Let  { S^,  i=0,l,2,...  } be  a renewal  process  with  inter- 
arrival cdf  F which  has  failure  rate  function  r(*)  , i.e.  F(x)  = 
x 

1 - exp(  - J r(s)ds  ) , x > 0 . Let  r be  a right-continuous  function 
0 1 

such  that  sup  r(s)  s r (t)  < 00  for  t 2 0 . Then  there  exists  a nonhomoa- 
CKs<t  1 

eneous  Poisson  process  { T^,  i=0,l,2,...  } with  cumulative  mean  function 
t 

/.^(t)  = r^fsjds  defined  on  the  same  probability  space  as  {s  , i=0 , 1 , 2 , . . . ' 

such  that  { Tp,  T*,  T*,  ...  ) 3 { SQ,  S.^  , S2>  ...  } almost  surely. 

Proof : As  is  frequently  done,  we  shall  not  explicitly  define  the  common 

probability  space  on  which  S and  T1  are  defined;  from  the  constructive 
definition  of  S from  T 


it  will  be  clear  how  an  appropriate  probability 
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could  be  defined. 

Consider  { T"*" , i=0,l,2,...  },  a homogeneous  Poisson  process  with 
1 

cumulative  mean  function  . We  shall  construct  a dependent  version 

of  { , i=0,l,2,...  },  the  desired  renewal  process,  by  "thinning"  T*  , 

i.e.  certain  points  will  be  removed  and  the  remaining  points  will 

constitute  { S^,  i=0,l,2,...  } . Define  SQ  = = 0 . If  = t^  , 


remove  ("thin")  the  point  with  probability  1-r  (t^) /r^  (t^)  , thus 

S1  = T1  probability  rft^/r^ft^)  conditional  on  = t7  . If 

is  thinned  and  = t^  < remove  the  point  with  probability 

l-r(t2)/r1(t2>  . Thus,  if  = ^ , T2  = t2  ' S1  = t2  with  Probability 


(1-r (t^) /r^ (t^) ) (r (t2) /r^ (t2) ) . Continue  inductively  until  the  first 


"unthinned"  point  of  T is  reached  and  define  this  point  to  be  . In 

other  words,  given  that  T = ( 0,  t,  , t_,  t,,  ...  ) define  an  integer- 

^ J 

valued  random  variable  : 

1 m / r<V  \ 

P ( M1  > m | T = (0,t1,t2,. ..)  ) =.Ri  [1  - r^(t j ; (2.3 


Then  S,  = T 
1 M 


We  now  consider  the  distribution  of  S,  . Let  N, (t)  = max{  n:  T1  < t } 

11  n 


P{  S > t } = E p{  S > t , N. (t)  = k } 
1 k=0  1 1 


p{  Si  > t | t1 tk  } f(t1,...,tk|k) 


0<t  <•  . .St,  St 
1 k 


dt.  . . .dt,  P{  N.  (t)  = k } 
1 k 1 


r \ k r,  (t.) 


= Z nil  - — ! k!  n -4 — - dt  . . .dt  X 

k=0  i=l  , r (t.)  1 i=l  A (t)  1 k 

i 11/  1 

0<t  <...<t  <t  \ / 

1 k 


exp(-A  (t)  ) A1  ^ 
k! 


Vfc 


n (r  (t.)-r(t.))  dt  ...dt  exp(  Ai(t)) 

i=l  1 i l 1 k 

k! 


= E ,nn  (r  (s) -r (s) ) ds  eXp(_Ai(t)) 
k=0  1"1  1 r-: 


-.S.  (A.  (t)  - A(t))k  exp(-Al(t)) 

k=o  i kT 


= exp ( A^(t)  - A(t))  exp ( -A^t)) 


exp(-A(t))  = exp ( - / r(s)ds 

0 
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The  third  equality  follows  from  Lemma  2 and  (2.3),  the  fourth  by  symmetry, 
the  rest  are  routine.  The  above  analysis  implies  that  the  distribution  of 
has  failure  rate  function  r(*).  Thus  has  the  desired  distribution 

Now  we  must  consider  the  succeeding  events  in  the  renewal  process  S . 
Condition  over  the  partition  { = T^  = s^  } , s^  > 0 ; the  random  index 


is  a stopping  time  ([3],  p.59),  therefore  the  strong  Markov  property 
([3],  p.131)  implies  that  the  distribution  of  { T^  . - T , i=0,l,2,...  } 

is  independent  of  { T'!’  , 0 < i < M,  } conditional  on  {t  = s,  }.  Given 

l 1 M l 

that  { T = S,  } , { T^  , . - T^"  , i=0,l,...  } is  a nonhomogeneous  Poisson 

1 M^+i 


process  with  cumulative  mean  function  A (t) 

J-  / J- 

Note  that  r (t+s  ) S sup  r(s)  £ sup  r(s) 

O^sSt+s^  0<s^t 

on  { T = s,  } , the  process  { T^  - T* 
1 M^+i 

same  hypotheses  which  enabled  us  to  construct 


= J r (s)ds  = / r (s+s  )ds  . 
sx  0 1 1 

for  all  t,  thus  conditional 

, i=0,l,2,...  } satisfies  the 

= T^  . Repeating  the 


constructive  procedure  gives  conditional  on  = s^  . Continuing 

inductively  allows  construction  of  the  entire  renewal  process  { , i= 

0,1,2,...  } . This  completes  the  proof  of  Theorem  1. 


At  this  junction,  some  remarks  are  in  order  concerning  "thinning": 

The  usual  connotation  of  "thinning”  refers  to  a more  random  form  of  deletinc 
points,  for  example,  for  0 < p < 1 , each  point  is  deleted  with  probability 
p independently  of  the  process  and  other  deletions.  It  is  well  known  that  a 
Poisson  process  which  is  "thinned"  in  this  sense  remains  a Poisson  process, 
see  Renyi  (1970),  p.254,  P.4.10.  Other  aspects  of  thinning  are  discussed 


i 


i 
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by  Rade  (1972)  and  Jagers  and  Lindvall  (1973)  . 

In  Theorem  1 if  sup  r(s)  < 00  then  letting  = sup  r(s)  the 

0<s<°°  0<s<°° 

hypotheses  of  theorem  are  satisfied  and  we  have  the  following  corollary. 

Corollary  1.  Let  { S^,  i=0,l,...  } be  a renewal  process  with  inter- 
arrival distribution  F which  has  failure  rate  function  r(').  If 

sup  r(t)  = r < 00  then  there  exists  a (homogeneous)  Poisson  process 
0<t  1 

{ T^,  i=0,l,...  } with  intensity  r defined  on  the  same  probability 
space  as  { S^,  i=0,l,...  } such  that  { T^,  , T^,...  } 3 { SQ,  , 

S^/...  1 almost  surely. 

Theorem  2 ■ Let  { , i=0,l,2,...  } be  a renewal  process  with  inter- 

arrival distribution  F which  has  failure  rate  function  r(’)  . Let 

r be  a right-continuous  function  such  that  inf  r(s)  5 r (t)  > 0 for 
0 s<t  0 

t 5 0 . Then  there  exists  a nonhomogeneous  Poisson  process  { T^  , i= 

0,1,2,...  } with  cumulative  mean  function  AQ(t)  = J0  r (s)ds  defined  on 
the  same  probability  space  as  { S^,  i=0,l,2,...  } such  that 

r 0 0 0 T r 1 

i T , T^ , T2,...  > C { S0,  S , S2,...  } almost  surely. 

Proof : Consider  { S^,  i=0,l,2,...  },  the  renewal  process  with  interarrival  I 

cdf  F which  has  failure  rate  function  r . We  shall  construct  a dependent 
r 0 , 

version  of  i T , i=0,l,2,...  } , the  desired  nonhomogeneous  Poisson  process, 
by  "thinnina"  f as  follows:  if  S.  = s.  and  S.-S  , = x. , delete  S.  with 

ill  j.-l  l l 

probability  1 - rQ (sj /r (xj , i=l,2 Let  T°  = SQ;  let  T^  equal  the 

0 

smallest  undeleted  S.,  i>l;  let  T_  equal  the  second  smallest,  etc.  We 

l 2 

r 0 , , 

use  Lemma  3 to  show  that  1 T^ , 1=0,1,...  ; is  a nonhomogeneous  Poisson 
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process  with  intensity  r (•)  • Let  N(t)  = sup{  n : S it},  and 

0 n n 

r 0 , 

N (t)  = sup{  n : T < t } 

0 n n 

Because  T°  is  a thinning  of  S,  the  O-fields  they  generate  are 
comparable:  a({  NQ(s),  s < t })  c o({  N(s),  s i t })  . Therefore,  in 

order  in  verify 


^ h"1  P{  NQ(t+h)  - Nq (t)  = 1 Nq(s),  s < t } = rQ(t) 


it  suffices  to  verify 


h ^ P { N (t+h)  - N (t)  =1  I N(s) , s i t } = r (t) 

IT^U  0 0 0 


By  the  strong  Markov  property  for  renewal  processes,  { N(s),  0 < s ^ S 


N(t) 


is  conditionally  independent  of  { N(s),  sN(tj  - s ^ t } given  S, 


N (t)  ' 


furthermore  since  S,_.  is  the  last  renewal  at  or  before  t , 

N It) 

0({  N (s) , S . . ^ s < t })  = 0({  S . , N(t)  })  ; also,  for  a renewal 

N(t)  N(t) 

process,  N(t+h)  - N(t)  is  independent  of  N(t)  given  sN(t)  • Putting 
all  of  this  together  gives 


P{  Nq (t+h)  - NQ(t)  = 1 | N(s),  s 5 t } 

OO 

= pt  N (t+h)  - N (t)  = 1 I N (t+h)  - N(t)  = k,  S„...  } * 
k=l  0 0 1 N(t) 

P{  N (t+h)  - N(t)  = k | SN  } . 


By  renewal  theory. 


P{  N (t+h)  - N(t)  = k ] Sv 


s*  > 5 


F(t+h-sfc)-F(t-sj> 
1 - F(t-sJ 


1 

i 
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equality  holds  for  k = 1,  strictly  inequality  for  k > 2 . Dividing 
(2.6)  by  h and  letting  h + 0 yields  rft-s*)  for  k = 1 and  0 
for  k > 1 . Thus  in  order  to  verify  (2.4)  using  (2.5)  , it  suffices  to 
show  that 


iiS pt  vtth>  - v11  - 1 1 N<t)  - i-  sN(t)  ■ •* 1 ■ 


V*1 

r(t-sj 


(2.7) 


But 


P{  Nq (t+h)  - NQ(t)  =1  | N (t+h)  - N(t)  =1,  S = s*  } 


t+h 

= ft  P(SN (t) +1  13  unthinned  I SN(t)+l  = S> 


dP{SN(t)+l  " S I 4 S SN (t) +1  " t+h'SN(t)  = S*} 


/, 


t+h  rft(s) 


t r(s-s#) 


s-s. 


r(s-s^)  exp(  - / r (u) du  ) 


ds 


t-s.  t+h-s. 

* . . . \ . # * 


exp(  - / r (u)  du  ) _ exp(  - / 

J 0 J 0 


r (u) du  ) 


t+h 

{ V8’ 


s_s* 

exp(  - / r (u) du  ) 


1 - exp(  - / 


t+h-s. 


t-s. 


r(u)du  ) 


ds 


t+h 

Ic  v> 


s-s* 

exp(-  / r(u)du  ) 
t-s 


1 - 1 + r(t-s*)  + O(h) 


ds 


-13- 


1 

h 


t+h 

K 


r0(5) 
r(t-s.) 


ds 


+ o (h) 


= Jl + o(h) 

This  verifies  (2.7)  and  consequently  (2.1)  of  Lemma  3.  Eq.  (2.2)  of  the 
Lemma  3 is  verified  similarly.  Thus  T^  is  a nonhomogeneous  Poisson 
process  with  intensity  rQ(  ) . 

Corollary  2.  Let  [ S^,  i=0,l,...  } be  a renewal  process  with  interarrival 

cdf  F which  has  failure  rate  function  r(  ) . If  inf  r(t)  = r > 0 

0<t<°° 

then  there  exists  a (homogeneous)  Poisson  process  { T*?  , i=0,l,...  } with 
intensity  r^  defined  on  the  same  probability  space  as  { , i=0,l,...  } 

, 0 0 0 i r i 

such  that  { TQ,  T^ , T — t c { SQ,  S^,  S 2’ ■ • • 7 almost  surely. 

We  shall  now  present  an  example  from  reliability  theory  which  indicates 
a possible  application  of  the  above  theorems:  bounds  on  the  expected 

discounted  cost  of  repair/replacement  of  a component  renewal  process.  Sup- 
pose the  life  time  of  a component  has  distribution  F with  failure  rate 
function  r(#)  . Let  S^,  S2>  ...  be  successive  failure  times.  Suppose 
if  a component  failures  at  time  t,  it  costs  c(t)  to  renew  it.  No 
assumptions  (except  integrability)  are  placed  on  c,  it  can  vary  daily 
or  seasonally,  etc.  Let  d be  a discount  factor.  The  discounted  cost  of 
maintaining  the  system  is 

OO 

. E..  exp(  -ds.)  c(S.)  . 
i=l  l l 
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0 1 . . 

If  T and  T are  nonhomogeneous  Poisson  processes  with  intensity 


functions  r^t*)  and  r (• ) respectively,  as  described  in  Theorems  1 


and  2 then 


.0,  ,0 


i|1  exp(-dTT)  c (tT)  exp(-dSi)  c(Si)  £2^  exp(-dT^)  c(T*' 


It  is  often  easier  to  evaluate  functionals  of  Poisson  processes  than 


renewal  processes.  In  the  above  case  it  is  much  easier  to  compute  moments: 


00  - dt  00  00 

/q  e ° (t)  r0  (t)  dt  < Egx  exp(-dS.)C(S.))  < fQ  e ”dt  c(t)r  (t)dt  . 


In  this  example  it  is  plausible  that  r(")  is  unknown  but  can  be 


assumed  i-o  be  bounded  by  r^  and  r^  thus  giving  the  above  bound  on 


expected  renewal  cost. 


3.  Comparisons  of  some  alternating  renewal  processes 


In  this  section  we  compare  some  processes  which  arise  when  failed 


components  are  not  instantaneously  repaired.  Consider  the  special  case 


of  a repairable  component  with  exponential  lifetime  distribution  F and 


NWU  (new  worse  than  used)  repair  time  distribution  G . The  component 


will  alternating  spend  random  times  functioning  and  under  repair.  Let  X 


be  the  performance  process  [2,2l]  : 


X (t)  = 


0 , component  functioning  at  t 


1 , component  under  repair  at  t 


If  Lp*  l2'...  and  Rj'  R2'*‘‘  are  t^le  successive  lifetimes  and  repair 


times,  let  T2i  = + R.)  and  TJi+1  = T2i  + Li+1  , i = 0,1,2 


J 
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X(t)  = 


0 T_,  < t < T 

2i  2i+l 


U 


T < t < T 

2i+l  2i+2 


i = 0,1,2,...  . Let  the  auxiliary  performance  process  [2] 


Z(t)  =' 


T < t < T 
2i  2i+l 


t-T  T < t < T 

2i+l  2i+l  2i+2  , 


i = 0,1,2,...  . This  process  describes  a system  which  is  functioning  at  time 

0 . Let  describe  a system  which  has  been  under  repair  at  time  0 for 

a duration  equal  to  u . Barlow  and  Proschan  (1976)  establish  a stochastic 

inequality  between  and  Z under  the  assumption  that  repair  times  are 

DFR  (decreasing  failure  rate) . We  shall  establish  an  almost-sure  inequality. 

Let  Rg  be  the  initial  repair  period  of  a system  under  repair  at  time 

0 . If  repair  commenced  at  -u  , let  TU^  = -u  and  = -u+RQ  > 0 . 

Define  TU  = + T.  , as  defined  above.  Then 

1 0 1 


zu(t) 


t-T 

2i+l 


u u 

T„.  < t < T , 
2i  2i+l 


u u 

T < t < T 

2i+l  2i+2 


i = -1,0,1,...  . Note  that  Z and  Z^  are  both  Markov  processes  with  the 
same  transition  law,  only  their  initial  state  differs. 


Theorem  3 The  processes  Z and  Z^  can  be  defined  on  a common  prob- 
ability space  such  that  Z(t)  S Z^(t)  , t a 0 , almost  surely. 

Proof : Let  R^  be  the  remaining  repair  time  in  a repair  period  which 
has  been  in  progress  u units  of  time.  Let  L^  be  lifetime  random 


1 
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I 


variables  and  R.  be  repair  time  random  variables,  i = 1,2,...  . Let 
Rfc  be  the  remaining  repair  time  in  a repair  period  which  has  been  in 
progress  t units  of  time.  Because  R is  NWU,  the  residual  repair 
time  R?  is  stochastically  greater  than  R^,  t>0  , i = 1,2,...  . By 
Lemma  1,  let  I = [o,l],  be  Borel  sets  and  be  Lebesgue  measure, 

define  R.  , Rfc  , t > 0 , on  (I.,  F p ) such  that  P(  R.  < r!"  ) =1  . 

li  ill  li 

Define  R^  on  (I„,  F_, P„)  and  define  L.  on  similar  probability  spaces 

0 000  l 

(J\  , G^,£L)  . Let  (ft,  F ,p)  be  the  product  space  of  all  the  above  spaces; 

extend  the  domain  of  the  above-define  random  variables  in  the  obvious  way. 

{ Z(t) , t > 0 } is  automatically  defined  on  (ft,  F ,P)  . We  define 

{ Z (t) , t S 0 } as  follows:  On  { R^  < L,  } , Z (t)  = u+t  for 
u 0 1 u 

0 < t < R^  , and  Z (t)  = Z(t)  for  t ^ rJJ  . Note  that  the  lack  of 
0 u 0 

memory  property  of  the  exponential  distribution  guarantees  that  the 


»■ 

i 


ii 

% 


first  lifetime  of  the  Z^-process  will  have  the  correct  distribution.  On 
{ Rq  > L1  } , condition  over  L^  = X,^  : if  L^  = X,^  and  R^  > X,  then 

at  time  X,  (when  the  Z-process  begins  its  first  repair  period)  the  Z - 
process  has  been  under  repair  for  a duration  u+X.^  , thus  R 1 

describes  the  repair  status  of  Z^  . On  the  set  { L^  = X,^  < R^  , 

rU+^i  < R + L } define  Z (t)  = u+t  for  0 5 t < X,,  + RU+^1  and 
112  u 11 

Z^(t)  = Z(t)  for  t 5 £ + R V1  . By  construction , the  Z-process  is 

repaired  before  the  Z^-process;  since  R^  1 < R + , the  Z-process 

will  be  functioning  when  the  repair  of  the  Z^-process  is  completed  and 

the  residual  lifetime  will  be  exponential.  This  procedure  of  defining  Z 


i 

t 


on  (ft,  F ,P)  can  be  continued  inductively.  The  next  step  is  to  condition 
over  { Li  = X^  , R1  = r1  , L2  = } and  { R^  > X^  } n { R^+X,l  > r = X„2  }. 
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Given  these  conditions,  the  Z^-process  has  been  under  repair  t^+u  time 

units  at  time  t^  = J^+r^+J.^  ' t^,us  R^+t3  t^le  distributi°n  °f  resid- 

r u 

ual  repair  time.  On  the  set  { > Rq  , = r ' L2  = ^2  ' 

R^+^l  5 r + £ , R^+t3  ^ R„+L_  ) define  Z (t)  = t+u  for  0 < t £ t_  + 

112223  u 3 

R^+t3  and  Z (t)  = Z(t)  for  t 5 t.+R^+t3  . It  should  be  clear  from  the 
2 u 3 2 

construction  that  Z^(t)  - Z(t)  with  probability  1 . 

Theorem  3 can  be  used  to  establish  almost  sure  comparisons  between  the 
operating  processes  of  coherent  systems  with  exponential  component  life 
times  and  NWU  component  repair  times  similar  to  the  stochastic  inequality 
established  by  Barlow  and  Proschan  (1976)  and  then  show  that  the  time 
until  first  system  failure  is  NBU. 


\ 


* 

! 
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